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Directed deterministic classical transport: symmetry breaking and beyond 
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We consider transport properties of a double 5-kicked system, in a regime where all the sym- 
metries (spatial and temporal) that could prevent directed transport are removed. We analytically 
investigate the (non trivial) behavior of the classical current and diffusion properties and show 
that the results are in good agreement with numerical computations. The role of dissipation for a 
meaningful classical ratchet behavior is also discussed. 

PACS numbers: 



I. INTRODUCTION 

Dynamical systems exhibit an extremely rich variety 
of behavior with regards to transport properties. For 
instance, it has been known since a number of years, 
that the chaotic nature of the dynamics may induce 
stochastic-like properties in a deterministic system, like 
normal diffusion, as it happens for a random walk |l| . A 
full understanding of how dynamics precisely determines 
the nature of transport in unbounded systems is however 
still not fully accomplished (for instance deterministic 
transport is observed also for non-chaotic, yet mixing, 
systems 0): moreover there is a wide set of systems 
for which such properties are quite subtle (for instance 
anomalous behavior may appear in systems with a mixed 
phase space as a result of long sticking times of chaotic 
trajectories in the vicinity or regular islands [3]). 

A property which has recently attracted much atten- 
tion is the "ratchet effect", namely the generation of 
transport with a preferred direction, in systems without 
a net driving force (or even against a small applied bias) 
4]. While stochastic ratchets are rather well understood, 
we are here interested in a purely deterministic setting, 
the starting point being in our case a family of area- 
preserving maps on a cylindrical (infinite) phase space. 
As remarked in [5j this generally requires a proper defini- 
tion of what is meant by "ratchet behavior" , as directed 
transport is easily achieved even for a free system, once 
the starting velocity is different from zero. A meaningful 
notion of ratchet behavior is that of getting a non-zero 
current for generic initial conditions without the action 
of a net force, provided fluctuations of the current are 
not so wide to overwhelm the effect. Such a behavior has 
been thoroughly investigated in [5[, for systems with a 
compact phase space: a sum rule has been established, 
allowing for a precise theoretical estimate of the ratchet 
current. 

We address however a completely different situation, 



where unbounded transport is not a priori ruled out by 
confining invariant structures, and the expected behavior 
is normal diffusion, with null average momentum. This 
choice is motivated both as it is theoretically challenging, 
in establishing examples of fully underdamped ratchets, 
and as modern cold atoms physics makes such systems 
good candidates for real experiments [f| 0, @ • Though 
the system under investigation may be easily quantized 
(see E3) we here address the classical setting only: 
our results will provide evidence that transport proper- 
ties, once we break up space-time symmetries arc quite 
non trivial, but a classical good candidate for ratchet be- 
havior requires an additional ingredient, which will be 
identified with dissipation in the present example. 



A. The model 

The paradigmatic example of transport in low dimen- 
sional Hamiltonian dynamics is provided by the Chirikov- 
Taylor standard map 



Pn+i = p n + ksm(O n ) 

On+l = n + Pn+1 



(1) 
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living in a cylindrical phase space (6,p) e S 1 x 1. Such 
a map does not only provide a case study to explore 
dynamical scenarios that arise as the nonlinear param- 
eter k is varied, but it also allows to consider novel fea- 
tures that appear upon quantization, most notably the 
so-called quantum dynamical localization Such a 

map exhibit symmetric transport properties: namely no 
average current ({pt — Pq) — 0) and linear growth of the 
variance (normal diffusion) in the regime k » 1 

As observed in the seminal paper [l3T |, searching for 
directed transport involves as a starting point a modifi- 
cation of the map ([T]) in such a way to break time-space 
reversal symmetry. This can be accomplished in a num- 
ber of ways [1, @, Q ; we choose to follow the suggestion in 
Q, as it is both theoretically simple and experimentally 
clean: namely unevenly spaced, phase shifted, kicks are 
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introduced, described by a potential of the form 

^,«(M= (2) 

cob (0) 

fc cos (0 - 0) £~ -nr- 1 + 0. 

The corresponding Hamiltonian is 

H M (p,6,t) = ±p 2 + V^(9,t). (3) 

p and 9 being conjugate variables. The corresponding 
discrete dynamics over a full period (corresponding to a 
pair of kicks) is written as 

' p n+1 = p n + ksm(9 n )+ 

k sm(9 n + £p n + £fcsin(0 n ) - 4>) 

< 

9 n +i = n +Pn + fcsin(0„)+ 

(1 - £)* sin(0„ + £p„ + £k sin(0 n ) - <j>). 

We immediately point out an important feature: the ex- 
tended map Q has an associated torus map, of size 2wM 
in momentum, only for £ = N/M, namely for commen- 
surate kicking times. 

Our primary goal will be to investigate transport prop- 
erties of the map ([3]), in particular the first two mo- 
ments {p t — pn) and ((p t — Po) 2 }, especially for a choice 
of parameters' values breaking space-time symmetries: 
in the present case this is accomplished by choosing 
£ ? {0, 1/2} and cj, + {0, n}. 

B. Transport 

First of all we observe that if we fix a pair £, <f> and 
vary the nonlinear parameter k, a standard map-like 
scenario appears (see fig. CO): namely for small k 
KAM invariant structures create a barrier to unbounded 
transport, but these are destroyed for larger nonlincar- 
ity, allowing in principle unbounded trajectories on the 
cylinder. 

Though transport may be anomalous if sticking regu- 
larity regions influence motion in the chaotic sea, the 
typical behavior is normal diffusion (linear growth of the 
variance), while symmetry breaking generally leads to a 
nonzero first moment. To be more precise we introduce 
the following notation: H^(t) will denote the first two 
moments at time t, obtained by averaging over a set 
Xo of initial conditions: in particular a natural choice 
is to consider initial sets A4 Po (po fixed, 9$ uniformly 
distributed on [0, 2%)). For a wide choice of parameter 
values we observe that, after a transient, ^m po (t) 

saturates to an asymptotic value Si , which has a 
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FIG. 1: Phase space surface of section. </> = n/2, £ = 0.2 and 
five different (increasing) values for the kick strength k. 200 
iterates; 300 initial conditions in (0, 2n) x (0, 2n/£). 

nontrivial dependence on po (see the full curve in fig. 
©)• 




FIG. 2: (p—po) after 50 couples of kicks, averaged over 100000 
initial conditions, versus po- k — 9, £ = 0.2, (f> — tt/2. Full 
line: numerical results; dashed line: analytic estimate. 

As pointed out in Q in this case we may obtain current 
reversal by tuning the phase shift between kicks; if we 
start from A^o, as a matter of fact, ITq changes sign if 
we go from (j> = <f> to (f> = —(f) (such a property is also 
consistent with our analytic estimates, as we will see in 
the next section), see the full line in fig. ([3|). 

So, by an appropriate choice of initial conditions we 
get an asymptotic momentum different from zero, and 
current reversal is easily obtained by tuning the phase 
shift between pairs of kicks; this can hardly be termed a 
ratchet behavior as the momentum distribution is broad, 
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FIG. 3: (p—po) after 50 couples of kicks, averaged over 100000 
initial conditions, versus <j>. k = 9, £ = 0.2, po = 0. Full line: 
numerical results; dashed line: analytic estimate. 



FIG. 5: {(p—po) 2 } — (p—po) 2 after 50 couples of kicks, averaged 
over 100000 initial conditions, versus po. k = 9, £ = 0.2, 
cj> — n/2. Full line: numerical results; dashed line: analytic 
estimate. 



with a diffusive spread (see fig. (j4|). 








FIG. 4: ((p - po) 2 ) - (p - Po) 2 , averaged over 100000 initial 
conditions, versus time (expressed in kicks couples), k = 9, 
$ = 0.2, = tt/2, po = 0. Full line: numerical results; dashed 
line: analytic estimate. 

Notice that also the diffusion constant exhibits depen- 
dence upon the starting set M. Po as shown by the full 
line in fig. ([5]). 

As observed in Q if the two kicks take place at close 
times peculiar effects may arise, in the sense that a cel- 
lular structure of the phase space emerges (see fig. ([S])). 

More precisely this feature originates by requiring that 
the product k ■ £ is small: to keep k£, small, we cannot 
lower k too much, otherwise the phase space is no more 
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FIG. 6: Phase space surface of section. (j> = 7r/2, f = 0.05 
and five different values for the kick strength k. 200 iterates; 
300 initial conditions in (0, 2?r) x (0,2tt/£). 



fully chaotic. The cells are of size 27r/£, separated by mo- 
menta p = ±(2m + l)7r/£ + where m — 0, 1, 2...; in- 
deed at these momentum values (called trapping momen- 
tum) Pn+i — Pn- They are not symmetric with respect 
to p — 0, reflecting the fact that, with broken temporal 
and spatial symmetry, we expect a net current; if <fi is 
equal to zero or n the cells are symmetric with respect to 
p = (no net current without breaking the spatial sym- 
metry). The darkest cell corresponds to the cell where 
initial conditions are located. 

Connected to this structure of the phase space, a 
typical trajectory spends a lot of time trapped in a cell 
before escaping onto another one, see figure (J7J). 
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FIG. 7: A typical trajectory (momentum versus time, ex- 
pressed in kicks couples); k = 13, £ = 0.05, = 7r/2. 



Transport is strongly dependent on initial conditions, 
in particular different behaviors emerge from initial 
conditions inside momentum cells or at cells boundaries, 
see for instance fig. |8]), where the asymptotic value of 
momentum is plotted. 




FIG. 8: (p) after 50 couples of kicks, averaged over 10000 
initial conditions, versus po- k = 13, £ = 0.05, <j> = it/2. 



We choose an initial condition near the trapping mo- 
mentum: pffl = (2m + l)7r + <j) + Sp and m = 0. We 
get: 

Pi ] = Pq ] + sin(0 o )+ 

sin(0 o + p^ + fc(«> sin(0 o ) - 4>) 

~ p ?) - (fc(«)) 2 cos(0 o ) sin(0 o ) - fc«%cos(0 o ). 
Exactly in the middle of the trapping region Sp = [l4| : 

or, if we turn back to the original variables: 

k 2 

Pi - Po ~ y£sin(20 o ) 
In a similar way for the angle variable: 
#1 = o +p a + ksm(e o ) + 

(1 - £)fcsin(0 o + p< fl + fe«) sin(0 o ) - 0) 

~ 0o +Po + £fcsin(0o)- 

(l-£)fccos(0 o )(<fy + fc^ sin(0 o )) 
that, when dp — 0, can be approximated to 

01 = 00 + P2 
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Here cell trapping makes the dependence upon initial 
data more transparent: unless we go too close to the 
cell boundary, averages don't depend on the exact initial 
condition but essentially only upon which cell we start on. 
As remarked in 0] , which we follow closely, near the cells 
borders the map can be approximated by the standard 
map (with a suitable kick strength) . Let us introduce the 
rescaled variable p& = £p and the effective kick strength 



FIG. 9: Local phase space, k — 4, £ = 0.05, <f> = 7r/2 (see 
text). 

In figure ([9]) we compare the local phase space around 
the momentum cell border in the double-kicked system 
(top picture) and in a sin(2x) single-kicked standard map 
with nonlinear parameter K = k 2 £/2 (bottom picture). 
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II. ESTIMATES FOR THE MOMENTS 

Analytic estimates for the first and the second moment 
of momentum distribution may be obtained by a stan- 
dard technique [l5j , consisting in a Fourier expansion of 
the (^-deterministic propagator, followed by inspection of 
leading contributions. If we start from an initial distri- 
bution P(p, 9, 0) we may write the average of the M-th 
moment as 



{(PN-P ) M ) 



(5) 



/ Q(O N , PN ,t N \O,p,O)P(O,p,O)(p N -p) M M N d0dp N dp 

where Q(6n,Pn, tN\9,p, 0) is the conditional probability 
of having 9n,Pn at time tjv given that at time we have 
8, p. In particular averages corresponding to M Po have 
an initial distribution 



P(9,p,0) = ±-6(p-p ). 
Q obeys the recursion property 

?N,PN,tN\6o,Po,0) = 



(6) 



(7) 

/ Q(6 N ,p N ,t N \6i,pi,ti)Q(6i,pi,ti\6 a ,p o ,0)d0idpi 

and may be expressed in terms of the deterministic <5 
propagators: 

Q(0i,Pi,U\di-i,pi-i,U-i) = (8) 

Et=-»«(Pi-ft-i + ^(M) 



(U - ti-i)(Pi-i - V"(0i_i)) + 2nki) 



where the sum over fcj occurs because 9 has the angular 
topology. 

Also {pn — po) is determined via map equations: 



itf-1 
1=0 



(9) 



The last ingredient we need is the Poisson summation 
formula, giving the Fourier transform of a 5-spectrum, 

hoo _^ -t-oc 

^2 ^(y + 27rn) = — ^2 exp[imy]. (10) 



n— — oo 



Note that the (^-function constraints S(pi — Pi—x + 
V'(9i-i)) take care of the p-integrals. To perform the 
integration we need the Jacobi-Anger expansion: 



Ji being Bessel functions of the first kind. 

The M-th order moment at time N, for our specific 
potential, may thus be rewritten as 

((PN-P0) M )= (11) 
mjv — — oo "' L^m\ — — oo Jo 2ir ' ' ' JO 2tt 

(/csin(0 o ) + fcsin((9i - <p) + ... + fcsin(0jv-i - 0)) M 

CXp [ i J2r=l m r(9 r - r -! - (t r - t r -\) 

(po + fcsin(6» ) + fcsin(6»i -</)) + ... + k sm(0 r _i*))) ] 

where r is a kick index, t r is the time and m r are integers; 
N is the last kick taken into consideration and we suppose 
it to be even. M is equal to 1 for the first moment and 
to 2 for the second one. 

(9 r _i* is equal to 9 r -\ if r is odd and to r _i — if r is 
even; similarly, t r — t r -i is equal to 1 — £ if r is even and 
to £ if r is odd. 

Since we are interested in the behavior of the map after 
the application of a couple of kicks, we define a new time 
variable: n = N/2. 

A. The current 

We want to estimate the current (p — po), i.e. equa- 
tion (fTTj) when M is 1. The so-called quasi linear result 
corresponds to setting all m,j — (fully random propa- 
gator) and is null in the present case. By taking only one 
to j ^ (note that it has to be ±1 in order not to have 
a vanishing integral) we get the first non zero contribu- 
tion to the average p. This contribution is referred to 
in the literature as one kick correlation contribution; it 
corresponds to the lowest harmonic contribution to the 
Fourier series for the <f> dependence of the current. 



(p-po)(n,M,O w ' 



(12) 



-fesnu>+(l-Opo)Ji[(l-£)fc] 
•A) 1(1 - K)k\ Wo[(1 _ e)fc] i + 



ksm((j> - ^Po)Ji[ik}\ 



-Mm 2 

(we have summed the terms coming from complex con- 
jugate (rrij — ±1) for every possible j-choice). 

The next frequency appears when we include the two 
kick correlations. This contribution is dominated by 
terms coming from evenly spaced kicks with coefficients 



of equal value, i.e. m,- = m 



j-2l 



±1. In the absence of 



any informations about the value of £fc, we take into con- 
sideration all contributions of this kind, summing over I 
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(for small £fc the dominant endowment is given by j and 
j — 21 as far as possible). We get 



(p-Po){n,k,<f>>0 



mj=±l,mj_2i=±l _ 



(13) 



fc(sin(20 + 2(1 - Opo) Ji[(1 - - 0*] 

(/i(Jo[(l-0*].Jb[2(l-0*D+ 

jo^ii-sjAJ i_j [2(i-e)fc]* ;+ 

fc(sin(-20 + 2^p )(Ji[^] 2 Ji[2^fc] 



(/ 2 (joKfc],j [2efc]) + 1 -^ffl 



L-Jo[2«fe] 

where /i and / 2 are given by 

, , „ a6 /l-(a 2 )"- 2 l-(fe 2 )"- 2 

Ji(a,o) = 

/a (a, 6) = 



a 2 - b 2 \ 1 - a 2 1 - b 2 

I I- (a 2 )"- 1 1 - (6 2 )™" 1 



6 2 



1-a 2 



1-& 2 



The next contributions would come from 



= ±1 and 



m j-(2i+i) = il or T 1 ! but their relevance is much 
smaller than the previous ones. 

Now we can compare these analytic results with the 
numerical data (figs. (0 [TUl fTTjl ). The agreement is 
satisfactory. 



v v vv 



FIG. 11: (p — po) after 50 couples of kicks, averaged over 
100000 initial conditions, versus £. k — 9, (f> — tt/2, po = 0. 
Full line: numerical results; dashed line: analytic estimate. 



momentum, and current reversal by a transformation on 
the phase shift (f>. If instead of a single valued initial 
momentum A4 Po we choose initial conditions distributed 
over the entire torus we have that the current goes to 
zero, both numerically and analytically. This can be 
done only for a rational value of £ (otherwise the torus 
unitary cell is not defined). For an irrational £ the 
current goes to zero if we calculate it over wider and 
wider intervals of initial conditions (see fig. (HH)). 



Aha, 



r 



W 1 



I 




FIG. 10: (p — po) after 50 couples of kicks, averaged over 
100000 initial conditions, versus k. £ = 0.2, <f> = 7r/2, po = 0. 
Full line: numerical results; dashed line: analytic estimate. 



FIG. 12: (p — po) after 50 couples of kicks, averaged over 
100000 initial conditions, calculated over wider and wider in- 
tervals of initial conditions, po is distributed in (—x/2,x/2) 
where x is the value in the abscissa-axis, k = 9, £ = 7r/10, 
<j> = 7r/2. Full line: numerical results; dashed line: analytic 
estimate. 



In particular we remark again that by properly choosing 
M. Po we may obtain non-zero asymptotic values of the 
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B. Diffusion 

In a similar way we can obtain analytic results for the 
behavior of the variance ({p — po) 2 ) — (p — Pa) 2 - 

The quasi linear approximation (i.e. taking all rrij = 0) 
gives Dql = fc 2 /4 (D is the diffusion coefficient: ((p — 
Po) 2 )-{ P -Po) 2 \ mj =o = 2D QL N) : 

Then we take only one rrij ^ 0; in this case we have a 
non zero integral by putting cither rrij = ±1 or rrij = ±2: 



rrij— ±1 



{(p-po) 2 )(n,k l( f>,O m > 
k 2 cos(^+(l-e)Po)((^o- Wo 



(14) 



l-Un 



w„ 2 



O T 2 l + ./o 2 -(2»-l)(.7 n 2 )"- 1 + (2n-3)(J 2 )" 

+fc 2 cos(0 - Cpo) ((Jo - Ji) 1 -^ 1 

l+(n-l)(J n 2 )"-»(J Q 2 ) 



- J 5- 



? (2J - 



(i-^) 2 



where the argument of the Bessel functions in the first 
part is [(1 — t;)k] and in the second part [£fc]. 



TUj — ±2 



((p~p ) 2 )(n,k^,O mj 



■^cos(20 + 2(1-Opo) 



(15) 



J a [2(l-0*]Jb[2(l-0*: 



l-(Jo[2(l-e)fc] 2 )"- 
1-Jo[2(l-Ofc] 2 



- £ cos(20 - 2^ ) J 2 [2£fc] ^gff" - 

The point now is to classify the contribution coming 
from a different choice of the rrij . With two rrij ^ the 
dominant contribution comes from evenly spaced kicks 
with coefficients of opposite value, i.e rrij — —rrij-21 = 
±1; the main endowment in this case is independent from 
the size of £k and is found by choosing terms as near as 
possible, namely the relevant contributions come from 
I = 1, 2, 3 . 



E i= l,2,3((P-Po) 2 )(«,fc,0,Cr j=±1 ' m - 2!=Tl = (16) 

(n - 2) + J [(l - Ok] 2 (n - 3) + J [(l - £)fc] 4 (n - 4)) 

-k 2 Mm 2 



( 



((n - 1) + JoimHn - 2) + J [^] 4 (n - 3)) . 

Summing up these contributions and comparing the 
result with numerics, we obtain a good agreement for 



the behavior of {(p — po) 2 ) — {p — Pa) 2 versus the initial 
momentum po, the kick strength k and the number 
of kicks couples n, see figs. (jU El [13]). About the 
dependence on k, notice that the contributions with 
some rrij 7^ only take care of the oscillations around 
an average value, given by the quasi linear result. 



A 




FIG. 13: ((p — po) 2 ) — (p — po) 2 after 50 couples of kicks, 
averaged over 100000 initial conditions, versus k. £ = 0.2, 
4> — 7r/2, po = 0. Full line: numerical results; dashed line: 
analytic estimate. 



III. DISSIPATION 

As we mentioned in the introduction the symmetry 
breaking paves the way for obtaining asymptotic non zero 
current, yet the linear growth in time of the variance 
broadens the momentum distribution, thus masking the 
asymmetry effect. 

A way to get a pristine ratchet behavior is to introduce 
dissipation y|, Q , while keeping the dynamics strictly de- 
terministic. 

To this end we modify the pair of maps in the following 
way: 



Pn 




jPn-i + fcsin(0 n _i) 

6 n -l + p n £, 

= jp n + k sin(6>„ - 4 

- e n + p n +i{i-0 



(17) 



(18) 



where < 7 < 1; if 7 = 1 we recover the hamiltonian 
system (no dissipation), while 7 = is the overdamped 
case. 
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We can also write equations (| 1 T|) and (fTH|) in the form 
of a single map including both kicks ((n — 1) — ► n): 

' p n+ i = 7 2 p„ + 7*;sin(6> n )+ 

k sin(6»„ + + £k sin(0„) - <j>) 

< 8 n+ i = 6 n + ^p n + (1 - S,h 2 Pn+ 

K + (l-07)*™(fl»)+ 

(1 - f )jfe sin(0 n + C7Pn + ^ sin(0 n ) - 0) 

In this case the presence of the time asymmetry £ is 
not required for the ratchet effect and it can be put equal 
to zero or to one half; in that case the ratchet effect is 
purely due to dissipation. On the contrary (f> is obviously 
still needed in order to break the space symmetry. 

It is known that very weak dissipation may lead 
to a quite complex organization of the dynamics (see 
[III), with many stable orbits with interwoven basins of 
attraction: deeper in the dissipative regime we typically 
observe, after a transient time, either an attracting 
periodic orbit, or a strange attractor (see fig. (jT4j) ) . 




e 



FIG. 14: Phase space surface of section: k = 9, 4> = 7r/2, 
£ = 0.2, 7 = 0.4. 

The most interesting features in this case are the 
following: i.) if space symmetry is broken the attractor 
is not generally symmetric in p direction, thus we expect 
a nonvanishing asymptotic current (see fig. (|15p ) 

ii.) In contrast to the conservative case here the finite 
size of the attractor prevents unbounded broadening of 
the distribution: after some transient time the variance 
saturates, see fig. (fT6)l . 

in.) Motion on the attractor seems ergodic, so the 
dependence on A4 Po here disappears, and (p n ) tends to 
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FIG. 15: Average current (p— po), averaged over 5 10 6 initial 
conditions, versus time, expressed in couples of kicks, po — 0, 
k = 9, £ = 0.2, cf> = tt/2, 7 = 0.4. 
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FIG. 16: ((p — po) 2 ) — (p — Po) 2 , averaged over 5 10 6 initial 
conditions, versus time, expressed in couples of kicks, po — 0, 
k = 9, ^ = 0.2, <t> = 7r/2, 7 = 0.4. 



an asymptotic value which seems independent on the 
initial probability density (see fig. (jTTJ) ) . iv.) Again by 
the (p — > — (j) we get current reversal (see fig. (HH)). 
In figure (|19p we plot the behavior of current versus 7. 



We are able to give an analytic estimate for the first 
moment also in the dissipative case, given some quite 
strict conditions on parameters: 7^1 (i.e. small dis- 
sipation) and £fe > 1. The last condition is related to 
the fact that when 7 is close to 1, the dissipative phase 
space is similar to the conservative one and the present 
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FIG. 17: (p), averaged over 100000 initial conditions, versus 
time, expressed in couples ol kicks for different values of po- 
Dot line: po = 0, dashed line: po = 20, full line: po = 100 
k = 9, £ = 0.2, (j> = tt/2, 7 = 0.4. 
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FIG. 18: (p), after 50 couples of kicks, averaged over 100000 
initial conditions, versus <f>. po = 0, k — 9, £ = 0.2, 7 = 0.4. 



approximation doesn't work if the phase space has a cel- 
lular structure. 

The general formula for the analytic first moment is 

(pn-Po)= (19) 

Z^m N = -oo ■■' {—im\=-oo JO 27r"'J0 2tt l^U /' 

fesin(0iv-i - 0) +7fcsin(0jv-2) + ••■ + 7 Ar ~ 1 fcsin(6» )) 

ex P i Y^=l m r [Or - 8r-l ~ (*r - £r-r)(??07 r + 

fcsin(6»*_ 1 )4- 7 fesin(6»*_ 2 ) + ... + 7 r - 1 /csin(6> )) 



FIG. 19: (p) after 500 couples of kicks, averaged over 10000 
initial conditions, versus 7. k — 9, £ = 0.2, = tt/2, po = 0. 



(notation as previously). 

The quasi linear result (i.e. setting all m,j = 0) is no 
more null, unless po — 0: 



<p-Po>Kfc,0,er j=o =Po(7 2 "-i)- 



(20) 



The next contribution is obtained by taking only one 
uij — ±1: 

(p-po)(n,fc,^O m ' =±1 = (21) 
-fcsin(0 + (1 - Opo) Ji[(1 - Ok]J [(l ~ 



^2n-3 l-((Jo[(l-Qfc]/7)) 2 r 

7 l-CJoKi-^feJ^)^ 



+fcsin( 



fOnl h rr/-l-,^n-2 l-((./oKfc]/7)) 2 )" 



Such an estimate reasonably reproduces numerical 
data for a fixed value of n (see fig. ([2H|) ). while getting 
meaningful asymptotic results is still an open problem. 



IV. CONCLUSION 

We have considered a two dimensional area preserving 
map, obtained by kicking a rotator twice, with the same 
strength, but with a phase shift. If such a shift, and 
the interval between successive kicks are chosen in such 
a way to break relevant symmetries, then the average 
momentum gain may differ from zero, but the broadness 
of momentum distribution hides the effect in the long 
time limit. The dependence of transport indices on the 
map parameters and on the initial distribution presents 
interesting features, which are studied both by numerical 
simulations and by analytic estimates. 
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FIG. 20: (p — po) after 500 couples of kicks, averaged over 
100000 initial conditions, versus 7. k = 9, £ = 0.2, <j> = ir/2, 
po = 0. Full line: numerical results; dashed line: analytic 
estimate. 



The most natural way to freeze (in a classical frame- 
work) the width of momentum distribution is to intro- 
duce dissipation (still keeping a strictly deterministic dy- 
namics): typically a strange attractor arises, and trans- 
port moments localize to a finite value which does not 
depend upon the choice of initial probability distribution; 
the phase shift can then be easily tuned to get current 
reversal. 
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